Ab initio study of compressed Ar (112)2: structural stability and anomalous melting 
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We study the structural stability and dynamical properties of Ar(H2)2 under pressure using first- 
principles and ab initio molecular dynamics techniques. At low temperatures, Ar(H2)2 is found to 
stabilize in the cubic C15 Laves structure (MgCu2) and not in the hexagonal C14 Laves structure 
(MgZn2) as it has been assumed previously. Based on enthalpy energy and phonon calculations, we 
propose a temperature-induced MgCu2 — >■ MgZn2 phase transition that may rationalize the existing 
discrepancies between the sets of Raman and infrared vibron measurements. Our AIMD simulations 
suggest that the melting line of Ar(H2)2 presents negative slope in the interval 60 < P < 110 GPa. 
We explain the origin of this intriguing physical phenomenon in terms of decoupling of the Ar and 
H2 degrees of freedom and effective thermal-like excitations arising from coexisting liquid H2 and 
solid Ar phases. 

PACS numbers: 31.15.A-, 64.70.K-, 81.30.-t, 64.70.dj 
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I. INTRODUCTION 

Solid hydrogen H2 is expected to become metallic at 
compressions higher than 400 GPa f]\. In fact, ex- 
perimental signatures of the long-sought insulator-to- 
metal phase transition remain elusive up to pressures 
of ~ 340 GPa 0. Accepted pressure-induced mecha- 
nisms by which the metallicity of hydrogen can be en- 
hanced involve atomization of H2 molecules and partial 
filling of electronic tr* molecular levels due to charge 
transfer from or band hybridization with other chemi- 
cal species Already in the earlier 90's Loubeyre 
et al., based on the disappearance of the Raman (R) vi- 
bron mode and the darkening of the material, claimed 
to observe metallization of the Ar(H2)2 compound when 
compressed in the diamond-anvil-cell (DAC) up to ^ 
200 GPa 0. The stable room-temperature (RT) phase 
structure of this compound was identified with the hexag- 
onal C14 Laves structure typified by the MgZn2 crys- 
tal (space group: PGS/mmc). Strikingly, posterior syn- 
chrotron infrared (IR) measurements did not show ev- 
idence of molecular bonding instabilities nor metallic 
Drude-like behavior up to at least 220 GPa @. Sub- 
sequently, Bernard et al. suggested that activation of 
H2 dissociation processes and corresponding develop- 
ment of metallization in Ar(H2)2 could occur via a solid- 
solid phase transition of the MgZn2 — >■ AIB2 (space 
group: P63/mmc) type at pressures already within the 
reach of DAC capabilities However, recent ab initio 
work done by Matsumoto et al. demonstrates that the 
onset of metallicity in the AIB2 structure commences at 
pressures significantly higher than in pure bulk H2 [Sj]- 

In view of the growing interest on hydrogen-rich van 
der Waals (vdW) compounds under pressure , partly 
motivated by the hydrogen-storage problem, and of the 
unresolved discrepancies described above, we have con- 
ducted a theoretical study on Ar(Il2)2 under extreme 
P — T conditions using first-principles density functional 
theory (DFT) calculations and ah initio molecular dy- 
namics simulations (AIMD). In this letter, we present re- 



sults showing that at low temperatures and pressures up 
to ^ 215 GPa the Ar(H2)2 crystal stabilizes in the cubic 
C15 Laves structure typified by the MgCu2 solid (space 
group: Fd3m). This structure has not been considered 
in previous works [tUsIIToI though its probable relevance 
to Ar(Il2)2 was pointed out recently ll']. On the light 
of first-principles enthalpy and phonon calculations, we 
propose a temperature-induced (pressure-induced) phase 
transition of the MgCu2 -> MgZn2 (MgZn2 MgCu2) 
type that may clarify the origin of the discrepancies be- 
tween the sets of R and IR data. Furthermore, in the 
high-P regime {P ~ 210 GPa) we find that a metallic 
hydrogen-rich liquid can be stabilized at temperatures 
of ~ 4000 K wherein H-H coordination features render 
molecular dissociation activity. By means of AIMD sim- 
ulations, we estimated an upper bound of the melting 
curve P(Tm) of Ar(H2)2 and found a negative dTm/dP 
slope spanning over the interval 60 < P < 110 GPa. 
Our simulations show that the lattice composed by H2 
molecules melts at temperatures significantly lower than 
the lattice of Ar atoms does, so leading to stable mix- 
tures of coexisting liquid H2 and solid Ar over wide P — T 
ranges. We propose an argument based on this atypical 
physical behavior to explain the cause of the estimated 
negative dTm/dP slope. 



II. OVERVIEW OF THE CALCULATIONS 

Our calculations were performed using the all-electron 
projector augmented wave method and generalized gradi- 
ent approximation of Wang and Perdew as implemented 
in the VASP code [l^- Dense Monkhorst-Pack special 
fc-point meshes fisj for sampling of the first Brillouin 
zone (IBZ) and a cutoff energy of 400 eV were employed 
to guarantee convergence of the total energy per parti- 
cle to within 0.5 meV. In particular, we used 8x8x8, 
8x8x4 and 12 x 12 x 12 fc-point grids for calculations 
on the perfect unit cell corresponding to the MgCu2, 
MgZn2 and AIB2 crystal structures, respectively. All 
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FIG. 1: Schematic low-T phase diagram of Ar(H2)2 under 
pressure. MgZn2 - MgCu2 and MgCu2 - AIB2 phase bound- 
aries are sketched according to the results and arguments pre- 
sented in the text. Thermodynamic states at which x-ray, R 
and IR vibron measurements were carried out are indicated: 
d x-ray = o, [3|] R = [l| IR = A and [l] IR = A. In- 
set: Enthalpy difference per particle of the MgCu2 (solid line) 
and AIB2 (dashed line) structures with respect to the MgZn2 
Laves phase as function of pressure at zero temperature. 
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FIG. 2: Calculated zero-temperature equation of state of 
Ar(H2)2 in the MgCu2 (solid line) and MgZn2 (dashed 
line) crystal Laves structures. Experimental data is from 
work [J./J^sef: calculated pressure-dependence of the H-H in- 
termolecular distance in the MgCu2 (solid line, solid dots) 
and MgZn2 (dashed line, solid triangles) crystal structures. 
Volumes and distances are in units of /particle and A, re- 
spectively. 



the considered crystal structures were relaxed using a 
conjugate-gradient algorithm and imposing the forces on 
the particles to be less than 0.005 eV/A. The phonon 
frequencies in our calculations were obtained using the 
small-displacement method f 14') 155 O'^^^ the unit cells (F- 
point phonon frequencies) and large supercells contain- 
ing 80 atoms. Ab initio molecular dynamics simulations 
were carried out in the canonical ensemble (TV, V, T) using 
bulk supercells of Ar(H2)2 containing 160 atoms (F-point 
sampling). At given pressure, the dynamical properties 
of the system were sampled at 400 K intervals from zero- 
temperature up to the melting curve of pure Ar. Temper- 
atures were maintained using Nose-Hoover thermostats. 
A typical AIMD simulation consisted of 3 ps of thermal- 
ization followed by 7 ps over which statistical averages 
were taken. It is worth noticing that we recently used a 
very similar computational approach to the one described 
here to study the energetic and structural properties of 
Ar(He)2 and Ne(He)2 crystals under pressure [ll|, and 
that very recent experiments |fl6] have confirmed the va- 
lidity of our predicted P{V) curves. 



III. RESULTS AND DISCUSSION 

A. Low-T results 

A series of candidate structures were considered in 
our ab initio enthalpy H{P) calculations (rutile, fluo- 
rite, MgNi2, etc.) however only the MgZn2, MgCu2 and 
AIB2 structures turned out to be energetically competi- 
tive so the following analysis concentrates on these ones. 



Ignoring quantum zero-point motion effects, we found 
that Ar(H2)2 is energetically more stable in the MgCu2 
structure than in either the MgZn2 or AIB2 structures 
at pressures P < 215 GPa (see Figure [1]). Also we pre- 
dicted that at Pt = 215(1) GPa the solid transitates from 
the MgCu2 to the AIB2 structure. The calculated zero- 
temperature equation of state of Ar(H2)2 in the MgCu2 
structure displays very good agreement with respect to 
the available experimental data H (see Figure [2]). For 
instance, at = 14.61 and 12.50 A^/particle we obtain 
a pressure of 8.83 and 17.14 GPa, respectively, to be com- 
pared with the experimental values 8.77 and 15.79 GPa. 
It must be noted that the estimated equation of state 
and pressure-dependence of the H-H intermolecular bond 
distance of the MgZn2 and MgCu2 phases appear to be 
indistinguishable within the numerical uncertainty in our 
calculations (namely, 1 GPa and 0.01 A-see Figure [2]-). 

Since hydrogen is a very light molecule, quantum zero- 
point motion corrections must be included in the calcula- 
tions. Customarily this is achieved using quasi- harmonic 
approaches that involve estimation of the vibrational 
phonon frequencies. In following this procedure, we 
found that Ar(H2)2 in the MgZn2 structure always ex- 
hibits imaginary F-phonon frequencies at pressures be- 
low ~ 300 GPa. This means that the MgZn2 structure 
is mechanically unstable, at least, at low temperatures. 
On the other hand, we found that the MgCu2 structure 
is perfectly stable at all the studied pressures. Our nu- 
merical tests showed that this result does not depend on 
the approximat ion of the exchange-correlation functional 
used [17| ■ In Figure [31 we show the phonon frequency 
spectra of Ar(H2)2 in the MgCu2 crystal structure calcu- 
lated along a /c-point path contained within the IBZ and 
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FIG. 3: Calculated phonon frequency spectra of Ar(H2)2 in 
the MgCu2 crystal structure at 1/ = 5.83 A^/particle and P ~ 
200 GPa (corresponding molecular vibron frequencies are not 
shown). The calculation was done on a supercell containing 
16 Ar and 64 H atoms and using a 4 x 4 x 4 fc-point grid for 
IBZ sampling. 



at pressure ~ 200 GPa . 

Since the x-ray resolved RT crystal structure of 
Ar(H2)2 is MgZn2 i, we carried out a series of 
AIMD simulations of this and the MgCu2 structure at 
T — 300 K in order to explore their stability. In 
fact, analysis of our simulations based on estimation of 
the averaged mean-squared displacement and position- 
correlation function [l^ showed that both MgCu2 and 
MgZn2 structures are mechanically stable at RT. On view 
of these results and the reported x-ray data, we concluded 
that a temperature-induced MgCu2 MgZn2 transition 
occurs at fixed P (or equivalently, a pressure-induced 
MgZn2 MgCu2 transition at fixed T - see Figure [1]-). 
This temperature-induced transformation can be under- 
stood in terms of entropy: Ar(H2)2 in the MgZn2 struc- 
ture is highly anharmonic so ionic entropy contributions 
to the total free-energy stabilize this structure over the 
MgCu2 phase with raising temperature. In fact, IR mea- 
surements displayed an splitting of the vibron mode in 
the interval 35 < T < 150 K This sppliting can be 

induced by the degenerancy removal of the infrarred Eiu 
mode due to the decrease of the crystal quality caused 
by precursor effects of the MgZn2 — >■ MgCu2 transition. 

The causes of the aforementioned R-IR experimental 
disagreements can be rationalized in the light of the fore- 
seen MgCu2 MgZn2 transition, as we explain in what 
follows. In Figure HI we enclose the calculated R and 
IR vibron frequencies of Ar (112)2 in the MgZn2, MgCu2 
and AIB2 structures as function of pressure. At high 
pressures, it is shown that the Aig (R) vibron line corre- 
sponding to the MgZn2 and MgCu2 structures are prac- 
tically identical, whereas they settle appreciably below 
the line estimated for the AlBo structure, and follow 
closely the low-T R results Q. These outcomes to- 
gether with the energy and phonon results already pre- 
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FIG. 4: Calculated R and IR vibron lines of Ar(H2)2 in the 
MgZn2, MgCu2 and AIB2 structures as function of pressure. 
Experimental data from j3|] (o) and jQ] (•) are shown. 



sented, led us to think that the sudden dissapearance 
of the R vibron mode observed at high-P might be re- 
lated to the MgCu2 MgZn2 transition unravelled in 
this Letter and not to pressure-induced metallic-like be- 
havior 3 or MgZn2 AIB2 phase transition as pre- 
viously suggested. In fact, accurate electronic density of 
states (DOS) analysis performed on the perfect MgZn2, 
MgCu2 and AIB2 crystal structures show no closure of 
the electronic band gap up to compressions of at least 
420 GPa. Furthermore, the Eiu (IR) vibron line that we 
calculated for the MgZn2 structure (see Figured]) agrees 
closely with Datchi's IR data obtained at high pressures. 
This accordance is coherent if one considers that the se- 
ries of IR experiments were realized near the edge of the 
MgGu2-MgZn2 phase boundary, as we sketch in Figured] 
Conclusive experiments confirming the validity of our 
statements might consist of new series of x-ray and IR 
measurements performed at temperatures below 100 K 
and high pressures; according to our predictions, the vi- 
bron signature will eventually disappear at the crossing 
with the MgCu2-MgZn2 phase boundary since IR vibron 
modes in the MgGu2 structure are inactive. 



B. High-T results 

An intriguing physical phenomenon has been recently 
predicted and subsequently observed in H2 under pres- 
sure. It consists in the appearance of a maximum 
peak on its melting line followed by a negative dTm/dP 



slope 20ll2l| . This phenomenon have been explained in 



terms of subtle changes on the intermolecular forces due 
to compression instead of more familiar arguments like 
promotion of valence electrons to higher energy orbitals 
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FIG. 5; Upper bound of the melting line of Ar(H2)2 under 
pressure (solid line). Melting states obtained in the AIMD 
simulations are indicated with • and corresponding error bars 
amount to 200 K. Thermodynamic states at which the lattice 
of hydrogen molecules is observed to melt in the simulations 
are represented by A (fitted to a dashed line). The mel ting 
line of pure Ar [24I (long-dashed and dotted line) and H2 [20| 
(dashed and dotted line) are shown for comparison. 
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FIG. 6: DOS of Ar(H2)2 calculated at P ~ 2i0 GPa (V = 
5.83 /particle) and different temperatures. Inset: Aver- 
aged radial H-H pair distribution function gu-H obtained for 
the liquid (solid line) and solid (dashed line) phases (distance 
is in units of A). Snapshots of liquid (L) and solid (S) config- 
urations as generated in the AIMD simulations; dissociated 
II2 molecules are observed in L. 



or occurrence of molecular dissociation processes. Moti- 
vated by these interesting findings on hydrogen, we in- 
vestigated the dynamical properties of Ar(H2)2 at high-P 
and high-T in the search of similar physical manifesta- 
tions and to provide further understanding of vdW com- 
pounds in general. There exist several well-established 
techniques by which one can determine the melting line 
of a material; these essentially base on solid-liquid phase 
coexistence simulations and/or Gibbs free-energy calcu- 
lations ^23"]. Application of these methods at quantum 
first-principles level of description, however, turns out to 
be computationally very intensive and laborious. A sim- 
ulation approach that has proved successful in reproduc- 
ing general melting trends in materials at affordable com- 
putational cost, including experimental negative melting 



nil. This 



slopes, is the 'heat-until-it-melts' method 
technique allows for estimation of a precise upper bound 
of the solid-liquid phase boundary of interest. It must 
be stressed that we did not pursue accurate calculation 
of the melting line of Ar(H2)2 but to identify possible 
anomalous effects on it. 

Results from our AIMD simulations are shown in Fig- 
ure [5] and can be summarized as follows: (i) the lattice 
of H2 molecules melts at temperatures significantly lower 
than the whole crystal does so leading to mixtures of liq- 
uid H2 and solid Ar over wide P — T ranges; (ii) the 
fusion of the H2 lattice occurs at temperatures very close 
to the melting line of pure hydrogen while Ar(H2)2 as 
a whole practically reproduces the melting behavior of 
pure Ar up to ~ 60 GPa; (iii) the value of the estimated 
dT„i/dP slope is negative within the pressure interval 
60 < P < 110 GPa whereas positive elsewhere. Re- 
sults (i) and (ii) can be interpreted in terms of similar 



arguments than recently disclosed in solid mixtures of Ne- 
He and Ar-He, namely: lattices composed of same-species 
particles effectively behave like not interacting one with 
another but following alike physical trends (bond dis- 
tance, compressibility, etc.) than found in their respec- 
tive pure system [11| . In fact, we recently conjectured 
the superior energy stability of the MgCu2 structure over 
MgZn2 in Ar(Il2)2, as rigorously demonstrated here, us- 
ing thi s ty pe of reasoning and crystal symmetry argu- 
ments [ill . Regarding result (iii) , this is in itself a man- 
ifestation of a very peculiar and intriguing physical phe- 
nomenon. In general, systems presenting negative melt- 
ing slope are characterized by open crystalline structure 
(water and graphite), surpassing promotion of valence 
electrons in the fluid phase (alkali metals) or continu- 
ous changing interparticle interactions (molecular hydro- 
gen). In the present case either open crystalline structure 
or electronic band promotion effects (energetically pro- 
hibitive) can be ruled out, so changes in the interparticle 
interactions holds as the likely cause. The question to be 
answered next then is: what is the exact nature of these 
changing interactions?, or more precisely, do they relate 
to dramatic changes in electronic structure? are the ionic 
degrees of freedom and effective decoupling of the Ar and 
H2 lattices crucial to them? In order to detect possible 
electronic phase transitions and/or molecular dissocia- 
tion processes, we performed meticulous averaged DOS 
and H-H radial correlation function (gH-n) analysis over 
the atomic configurations generated in the AIMD runs. 
Within the thermodynamic range 60 < P < 110 GPa 
and < T < 3500 K (where dT^/dP < 0), we found 
that Ar(H2)2 is always an insulator material, either solid 
or liquid, where H2 molecules remain stable. Therefore, 
dramatic changes in electronic structure can be discarded 
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as the originating mechanism behind the alternating sign 
of dT^jdP. Actuahy, on the contrary case, one would 
have expected the value of dTm/dP to be negative also 
beyond ^110 GPa since in principle there is not reason 
to think that the originating electronic effects go miss- 
ing at higher-P. Interestingly, at compressions above 
~ 210 GPa and temperatures high enough for the system 
to melt we did observe closure of the electronic energy 
band gap; according to gn-H analysis and visual recre- 
ation of the AIMD configurations, this electronic phase 
transition is a consequence of emergent H2 dissociation 
processes. In FigurelHl we enclose DOS and gH-H results 
that illustrate this insulator-to-metal phase transition. It 
is worth noticing that this result is consistent with shock- 
wave compression experiments [26^ and ah initio investi- 
gations |20|] performed on pure hydrogen. With regard to 
the predicted negative melting slope, ionic effects are the 
likely subjacent cause. Considering results (i) and (ii) 
above, one can envisage a plausible explanation of this 
phenomenon: since the lattice of H2 molecules melts at 
temperatures much lower than the lattice of Ar atoms 
does, collisions between diffusive H2 molecules and lo- 
calized Ar atoms act like effective thermal-like excita- 
tions on the last that ultimately provoke the global melt- 
ing of the system at temperatures below those of pure 
Ar. II2 itself presents negative melting slope beyond 
^ 60 GPa so this effect is echoed, and enhanced signifi- 
cantly, in the melting line of Ar(H2)2 (see Figure [5]). Be- 
yond ^110 GPa, the stability of the Ar lattice becomes 
further reinforced so the effect of H2 collisions there is 
just to deplete the slope of the global melting curve in 
comparison to that of pure Ar. In order to test the va- 
lidity of such hypothesis, we performed a series of AIMD 



simulations at P ~ 110 GPa {V = 7.50 A^/particle) in 
which we kept the lattice of H2 molecules frozen and left 
the lattice of Ar atoms to evolve. Under these condi- 
tions we found that fusion of the Ar structure occurs at 
temperatures ^ 900 K above that of pure Ar, hence the 
original melting mechanism proposed seems to be corrob- 
orated. Analogous melting behavior than described here 
can be expected in other rare gas-H2 compounds, like 
Xe (112)7, and maybe also in SiH4 (112)2 and CH4(H2). 



IV. CONCLUDING REMARKS 

To summarize, we have studied the behavior of 
Ar (112)2 under pressure at low and high temperatures us- 
ing computational first-principles techniques. As results, 
we have unravelled (i) temperature(pressure)-induced 
solid-solid phase transitions that may resolve the existing 
discrepancies between the sets of R and IR experimental 
data, and (ii) an anomalous melting phenomenon consist- 
ing of negative melting slope. The atypical melting line 
of Ar(H2)2 can be understood in terms of the decoupling 
of H2 and Ar ionic degrees of freedom and of coexistence 
of same-species liquid and solid phases. Metallization of 
liquid Ar(H2)2 is predicted at high P — T conditions. 
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